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ABSTRACT 


This  paper  extends  the  validity  of  exponential  tolerance  and  confl 
dence  limits,  under  certain  restrictions,  to  the  class  of  dlstrlbu 
tlons  with  monotone  failure  rate.  In  particular,  the  usual 
exponential  lover  tolerance  limit  la  shown  to  be  conservative  for 
the  increasing  failure  rate  class  of  distributions  in  the  range  of 
population  coverages  and  confidence  coefficients  of  practical  in¬ 
terest.  Conservative  confidence  limits  are  also  obtained  on  tall 
probabilities  and  moments. 


CONSERVATIVE  TOLLPuVJCE  /V.'O  CONFIDl.:CL  LINITS 


1,  Introduction 

A  fundamental  problem  in  statistical  reliability  theory  and  life  testir.>; 
is  to  obtain  lower  tolerance  Units  as  a  function  of  sample  data,  say 
X  »  ,  X^  t  ••••  ).  That  is.  If  X  denotes  the  time  to  failure  of  an 

item  with  distribution  F  ,  then  we  seek  a  function  K^v)  such  that 
P  {  1  -  F[L(X)  ]  ^  1  -  q  }  1  -  a  . 

We  call  1  -  q  the  population  coverage  for  the  interval  [L(X)  ,  and 

1  -  a  the  confidence  coefficient,  .'jnothe*'  important  problem  is  to  obtain 
a  function  N(X)  such  chat 

?  {  1  -  F(T)  ^  M(X)  }  i  -  c* 

for  a  specified  time  T  ^0  ,  Related  problems  are  those  of  obtaining  con¬ 
fidence  limits  on  moments  and  percentiles. 

Early  papers  in  life  testing  (e.g.  Epstein  and  Sobel  (1953)  )  derived 
confidence  limits  assuming  an  exponential  life  distribution.  Goodman  and 
Iladansky  (19o2)  examine  various  criteria  for  goodness  of  tolerance  intervals 
and  certain  optimum  properties  of  the  usual  exponential  tolerance  limits 
are  demonstrated.  Recently,  a  great  deal  of  effort  has  been  devoterl  to  ob¬ 
taining  various  confidence  lim.its  for  the  Weibull  distribution.  L'ubey  (L.  cJ 
obtains  asymptotic  confidence  limits  on  1  -  F(T)  and  the  failure  rate  for 
the  class  of  Weibull  distributions  with  non-decreasing  failun  rate.  e 
also  studies  the  properties  of  various  estimators  for  Weibi  11  parameters 


;nvbp'/  (1963)  ).  Johns  and  Lleberman  (1965)  present  a  method  for  obtaining 
i:.wer  confidence  limits  for  1  -  F(T)  when  F  is  the  Weibull  dtstri- 
rurton  with  both  scale  and  shape  parameters  unknown.  Unlike  Dubey,  they 
d:  not  require  that  the  Weibull  distribution  in  question  have  a  non-decreasf- 
ing  failure  rate.  These  confidence  limits  are  obtained  both  for  the  censored 
and  non -censored  cases  and  are  asymptotically  efficient.  Hanson  and  Koop- 
xuars  i'1964)  obtain  upper  tolerance  limits  for  the  class  of  dlst'-lbutions 
with  inc. ceasing  hazard  rate  and  lower  tolerance  limits  for  the  r.last=i  rf  dis* 
r-Lbutions  with  PF2  density,  f  (i.e.  log  f(x)  is  concave  v’here  finite). 
Hivever,  they  do  not  asstime  non-negative  random  variables. 

Assuming  that  the  sample  data  arises  from  a  distribution  with  monotone 
failure  rate  (either  non-decreasing  or  non-increasing  and  F(0~)  «  0  )  we 
obtain  conservative  confidence  limits  for  most  reliability  parameters  of 
interest.  These  confidence  limits  are,  in  part,  derived  from  the  exponential 
distribuciin.  Since  in  many  cases  these  are  optimum  confidence  limits  when 
t^e  failure  distribution  is  exponential  (Goodman  and  Madansky  (1962)  they 
arc.  in  this  sense,  best  possible  for  the  class  of  distributions  with  mono- 
rcne  failure  rate.  (See  Barlow  and  Proschan  (1965)  Chapter  2  and  Appendix 
2  fcr  a  discussion  of  distributions  with  monotone  failure  rate  and  a  test 
fir  its  validity.)  They  also  have  the  advantage  that  they  are  convenient 
(c  compute  and  are  not  based  on  a  strong,  non-verifiable,  parametric  assump- 
•  j'r..  Sin-.e  these  confidence  limits  are  derived  in  part  from  the  exponen¬ 
tial  disttibutlon  this  paper,  in  a  sense,  represents  a  new  justifiraLicn  fo> 
.-e  >f  exponential  confidence  limits  in  reliability  theory. 

?  Summary  and  Discussion  of  Results. 

let  X,  <X.,<...<X  <...  <X  denote  an  ordered  sample  from  a  itfe 

1—  2—  —  r—  —  n 
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distribution  F  .  We  shall  only  allow  the  possibility  of  censorship  on  tn 
right.  Our  methods  will  be  used  to  obtain  confidence  bounds  for  more  gen¬ 


eral  types  of  censorship  in  another  paper. 

We  say  that  a  distribution  F  is  IFR  (DFR)  if  and  only  if 

In  [1  -  F(x)  ]  is  concave  where  finite  (convex  on  [0,“]  ).  If  F  with  dt- - 

sity  f  is  IFR  (DFR)  then  the  failure  rate  f (t)  is  non-decreas *ng 

1  -  F(t) 


(non-increasing)  in  t.  Barlow  and  Proschan  (1964)  obtain  inequalities  for 
expected  values  of  statistics  based  on  the  exponential  assumption  when  in 
fact  the  true  distribution  has  a  monotone  hazard  rate. 

IFR  Results 

Let  .  r 


G 

r,n 


IX,  +  (n-r)  X 

,  i  r 


r 


and 


C,  (r) 
l-a,q 


-2r  In(l-q) 


if  Xi  (2r)  ^  -2n  In(l-q) 
i-a 


£ 

n 


if  Xi  (2r)  <  -2n  ln(l-a,) 

i-a  — 


where 


x:  (2r) 
1-a 


is  the  (l-a)-th  percentage  point  of  a  chi-square  distri¬ 


bution  with  2r  degrees  of  freedom. 

TllEORLM  1.  If  F  is  IFR,  F(0")  «  0  ,  sup  {  x  F(x)  ^  q  ' 

then 

(1)  P  {  1  -  F[C,  (r)  d  ]  >  1  -  q  }  >  1  -  a 


(2) 


'  C.  (r) 
-  l-a,q 


e  }  >  1 

r,n  - 


a  . 
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Mathematically  (1)  and  (2)  are  equivalent  statements.  When 

2 

>  -2n  ln(l  -  q)  •  the  lower  tolerance  limit  provided  by  (1)  is 

identical  with  the  exponential  tolerance  limit.  Amasingly  enough,  the  ex¬ 
ponential  lower  tolerance  limits  provide  conservative  tolerance  limits  for 
most  cases  of  practical  interest.  For  example,  if  l-a>l-e^'v  .633 
and  1  -  q  >  e  ,  then  the  Inequality  Xj_j,(2r)  >  -2n  ln(l  -  q)  holds. 

In  the  sense  of  being  "most  stable"  (see  Goodman  and  Madansky) 1962)  )  this 
Is  the  best  lower  tolerance  lisiit  for  the  exponential  distribution  and  hence 
a  "sharp"  conservative  tolerance  lisiit.  If  the  full  sample  is  known  this 
is  "best"  for  r  ■  n  . 


Let 


C*  (r) 
o,q 


(-lx  ln(l  -  q) 


X^2r) 


if  Xj^(2r)  <  -2  In(l-q) 


if  x‘(2r)  j:. -2  in(l-q) 


THEORLM  2.  If  F  is  IFR,  F(0")  -  0  ,  -  sup  {  x|f(x)  ^  q  } 

then 

(3)  P  {  F[C*  ^(r)  e,  „]  >  q  }  >  1  -  a 
a,q  r,n  —  — 


(4)  P  {  ;  <  C*  ^(r)  ^  }  >  1  -  a  . 

q  —  a,q  r,n  — 


In  this  case  the  exponential  upper  tolerance  limits  are  valid  when 

'f 

y‘'(2r)  ^  -2  ln(l  -  q)  .  Unfortunately  this  inequality  does  not  hold  for  all 
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values  of  r  (1  r  ^  n)  in  the  range  of  population  coverage  values,  q  , 
of  greatest  practical  Interest.  A  table  follows  which  gives  the  largest 
values  of  r  such  that 


X^(2r)  <  -2  ln(l  -  q)  . 


-  5  - 


Table  1 


Largest  values  of  r  such  that 
the  exponential  upper  tolerance 
limit  is  a  conservative  upper 
tolerance  limit  for  the  IFR  class. 

(i.e.  x^(2r)  <  -2  ln(l  -  q)  ) 

a  — 


6 


The  upper  tolerance  limit  given  in  (3)  is  a  significant  improvenient  over 
the  tolerance  limit  given  by  Hanson  and  Koopmans  (196A)  for  the  IFR  class. 
However  they  do  not  restrict  attention  to  non-negative  random  variables. 
Also  they  do  not  obtain  a  lower  tolerance  limit  for  the  IFR  class. 


THEOREM  3. 


If  F  is  IFR,  F(0")  -  0 


and  T  >  0  is  specified. 


then 

'5) 


P 


1  -  F(T)  >  6  (-  e  -  T)  exp 
—  n  r,n 


V 


t 


J  J 


>  1  -  a 


where 


«(x) 


if  X  ^  0 

if  X  <  0  . 


Johns  and  Lieberman  (1965)  study  the  problem  of  obtaining  lower  confidence 
limits  on  1  -  F(T)  for  the  Ueibull  distribution.  (3)  is  more  convenient 

r  * 

than  their  result.  However,  if  —  6  <  T  our  result  is  trivial.  If 

n  r  ,n 

r  “ 

—6  >  T  ,  then  it  is  identical  with  the  exponential  lower  confidence 

n  r,n 

limit.  In  reliability  applications  where  it  is  desired  to  establish  high 

reliability  the  mean,  hopefully,  will  far  exceed  T  and  therefore  it  seems 

quite  likely  that  6  will  also. 

r,n 


theorem  4.  If  F  is  IFR,  F(0")  -  0  and  0  -j,/”  x  dF(x)  , 
then 

(6)  P  {e<k  0  }>l-a 

—  oifT  r^n  • 
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where 


if  x^(2r)  >  2 

Q  — 


if  y^(2r)  <  2  . 

‘a 


Similar  confidence  limits  can  be  obtained  for  higher  order  moments. 

The  upper  confidence  limit  on  6  in  (6)  is  the  usual  exponential  con- 
2 

fidtnce  limit  when  <  2  .  Unfortunately  this  condition  is  not  satis¬ 

fied  for  values  of  r  greater  than  3  or  4  at  the  usual  significance  levels. 

In  acceptance  sampling  the  following  hypothesis  testing  problem  is  con¬ 
sidered  : 

H  ;  e  -  e 

o  o 

versus  H,  :  6  <  6  . 

1  o 

The  rejection  region  for  the  exponential  case  is  of  the  form: 

e  x^(2r) 

Reject  M  if  &  <  —  ~  — — 

o  r,n  —  2r 

2 

If  X  ^2r)  2  ,  then  by  (6)  this  test  is  also  a  size  a  test  for  the  IFR 

J 

<  e  ;  1  .  e  . 

p 


i  < 

r,n  - 


e  x^2r) 
o  y 


F  IFR 


0  > 
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DFR  Results 


As  we  might  expect,  if  a  useful  exponential  confidence  limit  exists 
for  a  problem  relative  to  IFR  distributions,  then  no  useful  exponential 
confidence  limit  exists  for  the  same  problem  relative  to  DFE  distributions 
and  conversely. 


THEOREM  5.  If  F  is  DFR.  F(0”)  «  0  ,  -  sup  {  x  )  F(x)  q  }  and 

Q^^r)  <  -2n  ln(l  -  q)  ,  then 


If 


-2n  ln(l  -  q)  we  can  only  make  the  trivial  statement 


P  {;  >0}>l-a. 

q  -  - 

For  most  cases  of  practical  interest  —  high  confidence  and  high  population 
coverage  —  (8)  is  not  a  useful  result. 


THEOREM  6.  If  F  is  DFR,  F(0”)  -  0  ,  -  sup  {  x  j  F(x)  q  }  and 

X^(2r)  >  -2  ln(l  -  q)  ,  then 


(10) 


-2r  ln(l  -  q) 
X^(2r) 


a 


-  9  - 


»■  11 )  p 


-2r  ln(l  -  q) 
xj(2r) 


3  V  ^  1  ~  “ 

r,nj  - 


The  iipj;er  ccnfidence  limit  is  trivial  when  x^(2r)  ^  -2  ln(l  -  q) 


Table  2 


Smallest  values  of  r  such  that 
the  exponential  upper  tolerance 
limit  is  a  conservative  upper  tol¬ 
erance  limit  for  the  DFR  class. 

(i.e.  x^(2r)  >  -2  ln(l  -  q)  ) 


I  -  01  *  .90 

1  -  a  -  .95 

1  -  a  -  ,99 

r 

r 

r 

1  i 

4 

4 

5 

.75 

4 

5 

6 

.80 

4 

5 

6 

» 8S 

5 

5 

7 

0 

.‘50 

5 

6 

7 

6 

7 

9 

.97 

7 

8 

10 

.98 

8 

9 

10 

8 

9 

11 

,  999 

11 

. .  -  - 

12 

15 

- J 

10  - 


THEOREM  7 


If  F  is  DFR,  F(0  )  ■  0  and  T  >  0  is  specified,  then 


(12)  P 


(- 


F(T)  >  6(T  -  re 


,  ^ 


where 


<S(x) 


if  X  >  0 


if  X  <  0 


as  before. 


THEOREM  8.  If  F  Is  DFR,  F(0")  -  0  ,  /  x  dF(x)  -  9  and 


X^(2r)  <  2n  ,  then 


(13)  P 


f  61  -  “ 


while  if  X*(2r)  ^  2n  ,  Chen 


(14)  F  if 

(13)  holds  for  significance  levels  of  practical  interest  when  r 


3.  Proofs  of  Theorems  in  Section  2. 

Let  Y  denote  a  random  variable  with  distribution  G  .  If  X 
continuous  distribution  F  ,  note  that  Y  ■  G  ^F(X)  has  distribut: 
We  will  repeatedly  use  the  following  lemma. 

Lemma.  If  G  ^F(x)  is  convex  non-decreasing  for  x  ^  0  ,  G  ^F(O') 
and  Y^  •  g”^F(Xj^)  ,  then 


has  a 
on  G 

«  0 


n 


r  ** 

r  ^ 

n 

-  a.  X, 

<  G 

I  a,  Y, 

1  .  ,  i  i 

Li=i  J 

(1) 

Li=i  J 

n 

wiien  a  ^0  and  I  a.  “1  (s,  or  a  *0,  i  =  l,2,...,n). 

i-1 

Tie  pru'if  is  obvious. 

In  what  follows  it  will  be  convenient  to  let 


G(x) 


/ 


1 


-X 

e 


V. 


0 


X  '  0 


X  <  0  . 


Proof  of  Theorem  1. 

Since  (1)  and  (2)  are  mathematically  equivalent  we  need  only  prove  (2). 
By  the  lemma  we  have 


I-  Y  +  (n  -  r)Y 

E  — 

i*l 


n 


>  F 


r  X.  +  (n  -  r)X, 


I 

i-1 


n 


since  G  ^F(x)  is  convex  when  F  is  IFR.  Now  choose  k,  so  that 

1-a 


P  <  (; 


r  Y,  +  (n  -  r)Y 

!  -i- 
1-1 


n 


<  k  [  -  1  -  a 
—  1-a 


2n 


i.e.  ln(l  -  k,  )  ■  Since  F  is  IFR  we  know  (Barlow  and  Prcschan 

i-a  - 


^bS) ,  p  27)  that 


0 


Kt;;  )  ^  b(t;c  ) 

q  -  q 


t/' 


1  -  (i-q) 


t  >  ; 


siiere  '  is  the  (unknown)  q-th  quantile 

q 


12  - 


Hence 


1  n 


>  F(  e  ;  ;  )  >  b(  ^  e  ;  ;  ) 
—  r,n  q  —  r,n  q 


and 


P  I  b  (  ^  0  ;  ;  )  <  k  }  >  1  -  a 

r,n  q  —  i-oi  — 


Since  b(t;  c,  )  is  non-increasing  in  ;  we  have 

q  q 


P  {  ;  b'^(^  e  ;  k- 

q  —  n  r,n  1-c 


)  }  >  1  -  a 


where  the  inverse  is  taken  with  respect  to  . 


Case 


k,  >  q  (i.e.  Xi  (2r)  >  -2n  ln(l  -  q)  ) .  From  the  follow- 

""  J.*Cl  ^  J. “ Cl 


ing  figure 


-2r  ln(l  -  q)  0 


^q  - 


r .n  . 


Case  2.  k  <  q  (i.e.  x,  (2r)  <  -2n  ln(l  -  q)  ) .  In  this  case 
X  *  vJ  X  ”  S.A 


b(  77  c  ;  s  )  k,  if  and  only  if  ;  ^  x 

"  r,n  q  —  1-j  q  ^  r,n 

In  either  case  we  have 


Pi-  -  C,  (r)  •• 

q  -  I-  q  r,n  - 


-  >3 


I 


if  Xi.„(2r)  >  -2n  ln(l  -  q) 

if  Xi.^(2r)  <  -2n  ln(l  -  q) 

Pro  n  of  Theorem  2.  Again  we  need  only  prove  statement  (4) .  We  use  the 
tcLlovlng  inequality  which  follows  from  the  IFR  assumption  and  the  lemma: 

r  r 

We  'hocse  k  so  that 
a 

P  {  G  I  Z  Y^  +  (n  -  r)  Y^  ]  L  \  1  « 

2(2  V 

je.  In(l-k)*  .  From  Barlow  and  Proschan  (1965)  p.  27  we  have 

‘  s  ■ 

g  . 

2 

Case  I.  k^  >  q  (i.e.  Xg,(2r)  ^ -2  ln(l  -  q)  ) .  From  the  following 

A  A 

figj'-e  we  see  that  B(r  0  ;  c  )  ^  k  if  and  only  if  C  <  r  0 

r^n  q  g  q  r^n 

“  14  * 


B(r  6  ;  c„)  >  k  If  and  only  If  ;;  '  -2r  in(i  -  9) 

r.n  d  -  a  q-  2 

a 


■II 


Cases  1  and  2  together  establish  statement  (4) 
Proof  of  Theorem  3.  Again  we  use  the  inequality 


[Z  Y.  +  (n  -  r)  Y  I  I  I  X.  +  (n  -  r)  X  / 

- -J-  - ':! 


and  choose  k,  so  that 
1-a 


(  f  Z  Y  (n  -  r)  Y  1  ( 


i.e.  ln(l  -  k  )  -  .  Let  p  -  F(T)  .  We  again  use  the 

2n 


sharp  bound 


Co 


F(t;  p)  ^  b(t;  p)  • 


t  <  T 


1  -  (1-p)^^^  t  il  T  - 


Then 


’■  1  “  »  *t.n*  P'  -  “l-o  I  -  ^  ■ 

Since  b(t;  p)  is  increasing  in  p  ,  we  have 

P  I  1  -  F(T)  >  1  -  b-1  (i  J  1  1  - 


-  15  - 


(5)  follcws  when  we  recall  that 


0  wiien 


^  ^.n‘ 


-  e  ,  I 

n  i-.x! 


P.cQof  of  Theorem  4.  To  show  (6)  use  the  sharp  bound 


1  -  e 


-t/e 


F(t;  b)  ^  B(t;  6)  * 

.Barlow  and  Pioschan  (1965),  p.  27) 
together  with 


t  ■■  G 


c  ■-  e 


G  i  Y,  +  (n  -  r)  Y  ]  <  F(r  6  ;  6)  <  B(r  0  ;  0) 

i  r  —  r,n  —  r,n 

to  estahlish 

P  {  B(r  0  ;  6)  >  k  }  '  1  -  a 


r.n 


—  a  — 


where  inti  -  k  )  *  -y  (2r)  as  before, 
a  a 


-1  2 

Case  1 ,  k  >  1  -  e  (i.e.  x  (2r)  >  2).  From  the  following  figure 


wfc  see  that  B(r  9  ;  9)  >  k  if  and  only  if  0  <  r  0 


r.n'  -  a 


r.n 


-  16  - 


Case  2 . 


-1  2 

k  <  1  -  e  (i.e.  x  (2r)  <  2).  Also  we  see  that 


B(r6  ;G)>k  if  and  only  if  -  <  r,n 

r,n’  a  — — • — 

.ience  the  result 

X  (2r) 

a 

To  show  (7).  Use  the  sharp  bound 

(“ 

t  < 

F(t;  ^  b(t;  =  J 

1  -  e 

t  - 

where  w  depends  on  t  and  satifies 
.t  -wx  , 

(16)  ^  e  dx  =  V  , 

n 

o 

(see  Barlow  and  Proschan  (19o5)  p.  28), 
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is  a 


which  establishes  (7). 

We  omit  proofs  of  the  DFR  results  since  they  are  a  straightforward 
application  of  the  same  techniques  applied  to  bounds  on  DFR  distributions. 
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